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We investigate the possible superconducting pairing symmetries which would account for the 
latest experimental results on the Bechgaard salts. Using a renormalization group (R.G) technique 
we calculate the couplings of the singlet and triplet interactions in the quasi ID Hamiltonian. In 
the presence of interchain interactions, the singlet and triplet couplings are of the same order of 
magnitude. The renormalized couplings are then used in an RPA calculation to evaluate the critical 
fields of superconductivity for the two dominant order parameter symmetries in the R.G flow, i.e, 
the d x 2_ y 2 and f y symmetries. It is shown that, for the standard values of the anisotropy ratio 
of the Fermi surface, the critical fields of the singlet symmetry are strongly reduced. A reentrant 
superconductivity is however still present for the triplet state. 
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Introduction Since the discovery twenty eight years 
ago of superconductivity in the organic compounds of the 
(TMTSF) 2 X, family ( X= PF 6 , ClO A , ) Q, Q, these 
compounds have been intensively studied both theoreti- 
cally and experimentally. Many original phenomena has 
been discovered and studied. However, at the time be- 
ing even the symmetry of the superconducting state is 
not yet known with certainty. The simplest question of 
which pairing symmetry characterizes the order parame- 
ter is still a matter of debate. Many apparently contra- 
dicting experimental results have attempted to address 
this question. Many have suggested the spin-pairing of 
electrons to be of the singlet type, while others suggested 
the spin-triplet. Our purpose, in this paper, is to study 
which theoretical model for the superconducting state is 
able to account for all the experimental data. 

At zero magnetic field, the effect of non magnetic im- 
purities 3 has been shown to reduce the critical temper- 
ature, in agreement with an unconventional symmetry, 
and with the presence of nodes of the gap function on 
the Fermi surface [4|]. 

Measurements of the critical field in the a and b' 
[E @, 0i 9 directions performed on both PF 6 and CIO4 
compounds have reported a superconducting state sur- 
viving up to fields as high as 9 and 5 T respectively. In 
a and b' directions, the maximum critical fields reached, 
are two times higher than the Pauli paramagnetic lim- 
iting field. At high magnetic fields ~ 4T, an increase 
in the critical temperature in increasing magnetic field 
has been previously observed, which was ascribed to a 
reentrant superconducting phase 0|- 

Lee et al. measured the electron spin contribution to 
the total Se Knight-shift in the PF§ compound for 
fields parallel to the a and b' directions respectively Q. 
In both cases it was found that the resonant frequency 
shift remains unchanged when the system passes through 
the superconducting transition, whereas it is expected 
to decrease and ultimately vanish if the pairing were of 
the singlet type. Later, Shinagawa et al. [l(J performed 
analogous measurements on CIO4 compound for a field 



strength of 1.38 T, and drew the same conclusions. How- 
ever, very recently, new Knight-shift measurements per- 
formed by Shinagawa et al. on the CIO4 at lower field 
values (0.9 T) in the a and b' directions, showed a de- 
creasing spin susceptibility upon cooling through the su- 
perconducting transition, while it remained unchanged 
for a field of 4 T. 

On the theoretical level, there have been three main 
ideas proposed to interpret the experimental results. The 
first idea was that the pairing symmetry of the order 
parameter is singlet at all magnetic fields. The high 
field re gim e is only marked by a first order transition 
[TTI . [l2l Il3| from a homogeneous superconducting state 
to a LOFF state [M [IE 13 G3 where the Cooper pair 
has a finite total momentum in order to compensate for 
the Zeeman effect. In the second interpretation, the 
superconducting state is assumed to be in the triplet 
[II GJ, [H Equal Spin Pairing (ESP) state at all mag- 
netic fields. Using this symmetry, with an anisotropic 
d-vector order parameter, Lebed et al. [l8| could explain 
the high field Knight-shift results. It has also been sug- 
gested [2H . l22l l23l |24| that in an increasing magnetic 
field the symmetry of the order parameter could change 
from singlet to triplet so that the superconductivity is no 
longer paramagnetically limited. In our previous work 
[23 | . we could calculate the transition line between a d- 
singlet and a p-triplet superconducting state using the 
free energy criteria. There, we considered the Zeeman 
effect as the only limiting effect on superconductivity, 
based on the fact that at the transition field the con- 
finement of the electron orbital motion to the conduc- 
tion layers reduces considerably the orbital effect of the 
magnetic field on superconductivity. This field induced 
dimensional crossover has been discovered by Lebed in 
his pioneering work of ref. [Ifl, which is essential to the 
physics discussed here. 

In the present paper we study a d x i_ y i singlet and an 
f y triplet state at the vicinity of the critical field, when 
the field is parallel to the b' direction. First, we will start 
our calculation by evaluating the renormalized triplet 
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and singlet couplings using a renormalization group ap- 
proach. Then, we will use the renormalizcd couplings in 
the calculation of the critical fields, for both symmetries 
of the order parameter. 

Theoretical model In order to explain the ID behav- 
ior of the Bechgaard salts at high temperatures, it is 
well known that the coupling between Cooper channel 
(electron-electron pairing) and Peierls one (electron-hole 
pairing) has to be taken into account. This is possible by 
means of a renormalization group (RG) approach [26J, 
used to calculate the effect of the high T ID fluctua- 
tions on the different couplings. When temperature is 
decreased, the system undergoes a dimensional cross-over 
at T = T* ~ 300A" where its 3D character becomes pre- 
dominant [26J] . The Cooper and Peierls channels are then 
decoupled, and the RG flow may be turned into a stan- 
dard RPA calculation. When a magnetic field is present, 
a three cut-off RG approach [2?], HII is more appropri- 
ate. In the following we use this RG method in order to 
determine the renormalized singlet and triplet couplings. 

At high temperature, we use the standard g-ology 
model, with hamiltonian H — i/kin + ^int , where kinetic 
and bare interaction parts write 
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interchain scattering, which operates as a modulation of 
the intrachain ones. 

The three cut-off RG method involves three charac- 
terictic energies, fisH -C fc^T* <C A (A is the energy 
band-width). For T > T* , the system is effectively ID, 
and the energy scale of the magnetic field as well as that 
of t c are irrelevant. 

We use the One Particle Irreductible (OPI) scheme 
[2!j[3(],[3l|,[32j to recalculate the scattering couplings in a 
one-loop expansion. The main idea of this RG method is 
to sum up high energy momenta in the action in order to 
get an effective action in terms of low energy momenta. 
This is done by reducing the energy scale A = A e~ e , 
during the flow (£ is the flow parameter and the flow be- 
gins with A = Ao). For each A, we get effective couplings 
<7i(A) (i = s,t) while gi(Ag) are the bare couplings. Using 
the OPI scheme, we write the RG equations relating the 
^■(A) to the & (A). The lengthy RG equations will not 
be listed here, for details see Ref. [29|, [3(| . In our numer- 
ical calculation of the renormalized couplings, there are 
three independant adjustable parameters, i&, the bare 
values of gi (where we chose bare couplings g s ~ 3gt 
which does not affect the final point in the flow) and the 
inter-chain backward scattering coupling We have 
used the following values g s /(nvf) = 1.04, g t /(nvf) ~ 0.4 
and gtb/i^Vf) = 0.038 which are compatible with known 
experimental data. At A = A* ~ tb, which fairly corre- 
sponds to T ~ T* , we get the renormalized couplings : 
g s * = 0.361 and g t * = 0.312. For T < T*, the couplings 
between Cooper and Peierls channels become negligible, 
so that the RG equations are equivalent to those obtained 
in RPA. Because of the decoupling of the RG equations 
for T < T*, the values obtained for the renormalized 
couplings at A* remain unchanged for T < T*. 

If we continue the RG flow up to A c , where some cou- 
plings begin to diverge, we get the following singlet and 
triplet response functions (see Fig. [1]) , where one can ob- 
serve the competition between triplet /-pairing suscepti- 
bilities and singlet d-pairing ones. In fact, by increasing 
the CDW fluctuations, the inter-chain backward scatter- 
ing enhanced the /-triplet superconductivity, resulting in 
renormalized triplet couplings of the same order of mag- 
nitude as those of the singlet ones. 

Renormalized couplings, g s * and <? t *, will be used in 
the following low temperature effective mean field hamil- 
tonian. 



where cs = cos((k y + k' ).b). 

The kinetic Hamiltonian is linearized around the Fermi 
points with a Fermi velocity Vf. a, b and c are the lattice 
parameters in the x, y and z directions, where for sim- 
plicity we assumed an orthorombic structure, a = (+/— ) 
stands for left and right moving electrons respectively, a 
is the electron spin. 

In the interaction Hamiltonian, the first three terms 
are related to intrachain interactions, while the last three 
ones are related to interchain interactions. In these inter- 
chain terms, g s and g t are replaces by gbb, the backward 
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where singlet and triplet order parameters are 



A fl (k,q) = A u (k,q) 



(3) 



3 




j,(px) 



(py) 



X 1 / 
z (tx) 



Figure 1: Flow of the superconducting susceptibilities \^ 
(singlet of symmetry d), Xt^ (triplet of symmetry p x ), 

X( Py ' (triplet of symmetry p y ), Xt^ (triplet of symmetry 

( t ) 

f x ) and Xt v (triplet of symmetry f y ), with bare parame- 
ters g s /(irv f ) = -1.04, gt/(nvf) = -0.4, g bb /(wv f ) = 0.038 
and 2t b /A = 0.133. 
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k is related to the internal degrees of freedom of the 
pairing state and gives the symmetry of the wave func- 
tion, while for the singlet order parameter q is related 
to the motion of the center of mass of the Cooper pair 
and gives the spatial variations of the gap function, q is 
therefore the LOFF wave vector. 

The d x 2_ y 2 and f y symmetries, and the corresponding 
interaction channels, are given by 



A s (k,q) = A s (q) cos(fc ± ) (5) 
g s *(k,k') = g s * cos(k±) cos(k' ± ) 



A t (k) = A t frin(2*j.) 
flf a *(k,k') = g s * sin(2fc ± ) sin(2k' ± ) 



(6) 

(7) 
(8) 



k X f a {k y 



V X fjky 



= kf 

<7/ifcg 

v f 



^y- cos(bk y ) + cos 2 (bk y ) 



cos(6 ky)] 



2at a sin [ak x fl7 (k y )\ 



kf a (ky) is the x component's module of the Fermi wave 
vector of the spin a electrons, and v^ a .(ky) is the corre- 
sponding Fermi velocity. (3 = v2 tb/t a . /!& is the electron 
magnetic moment, a, b and c are the lattice parameters. 

Unlike the above RG calculation, in the low tempera- 
ture regime, i.e below the critical temperature for super- 
conductivity T c (0) ~ IK, the magnetic field is relevant 
and should be considered in the mean field hamiltonian 
Eq.([2]). Its Zeeman effect is already included through 
the spin dependance of the energy spectrum. The or- 
bital effect will be fully taken into account through a 
Peierls substitution, which is justified in our case. More- 
over, Eq.([n]) takes into account the real quasi-ld form of 
the Fermi surface as was first noticed by Lebed [Hj], in 
contrast to previous calculations where a flat ID Fermi 
surface was used [ill, E3, E HI HI . 

In ref. [13] we have studied the hamiltonian Eq.Q 
with the Zeeman splitting as the only effect of the mag- 
netic field. There, we neglected the transverse transfer 
integral t c compared to t a and tb, and therefore neglected 
the orbital effect of the magnetic field. By calculating the 
gap equations, we could study the d x 2_ y i and p y symme- 
tries of the superconducting order parameter, p, d and / 
symmetries are believed to be the most probable candi- 
dates for the symmetry of the superconducting state in 
the Bechgaard salts [l|, H, [13, HH . The ga p eq uations for 
the / order parameter and for the p one [33| are identi- 
cal. Using the free energy criteria, we have been able to 
calculate a first order phase transition line between the 
d-singlet and the p-triplet superconducting states. How- 
ever, in our calculation of the phase transition we did not 
take into account the possible occurence of a LOFF state 
which should succeed to the homogeneous singlet state 
in increasing magnetic fields. Moreover, from the above 
RG calculations it appears that /-wave symmetry is to 
be considered rather than the p one. 

In this paper, we will establish the gap equations for a 
d x 2_ y 2 and a f y order parameters in a magnetic field in 
the b' direction. 

The electron Green's function in the b' direction is 
given in the gauge A = (0, 0, —Hx) [H, [3?J by 



We rewrite the linearized energy dispersion Eq.JT]) as 
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where Fermi velocity and Fermi wave number now de- 
pend on the transverse position on the Fermi surface k y 
as 
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and A c = it c /(vfG) is the semi classical amplitude of the 
electron motion in the intcrplane direction c [251 ] . 

Using the equation of motion for the above Green's 
function we can write the linearized gap equations for 
(T - T C (H)) < T C (H) as follows 
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for the triplet order parameter. 

d c is a cut-off distance. = (2^ b H/vf)0 . g Stt * 
are the dimensionless interaction constants determined 
in the RG calculation. The singlet superconducting 
gap function was assumed to have the form A q (k, r) = 
A q (r) cos(fcjJ. And we anticipated the fact that the 
maximum critical field corresponds to the LOFF wave 
vector perpendicular to the field direction, so that 
A q (r) = &.q x {x) = A cosmic). It should be noted that 
the /y-gap equation (TT2"]) is identical to the p y -gap equa- 
tion calculated previously [33|. The only difference re- 
sides in the values of the p and / couplings. 

Results and discussion The integral in the /-wave gap 
equation (fT2")l diverges for T — > 0, indicating that there 
is no finite value for the critical field. This is at the 
origin of the reentrance phase of superconductivity at 
high magnetic fields for the triplet case. On the contrary, 
the integral in the <i-wave gap equation (fTTj) converges 
at T = 0, and therefore the zero temperature critical 
field for the singlet phase exists and is finite, this is the 
paramagnetic limit of the LOFF state, i.e H^ OFF . 

Let us first revisit the case without orbital effect, where 
t c = 0. Taking the limit T -> in Eq. ^ we find 
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Figure 2 
of the d a 

dashed line shows the critical field for the homogeneous state, 
while the solid line shows the LOFF-state critical field after 
maximization with respect to the LOFF wave vector. T cs ~ 
I.IK is the critical temperature corresponding to the singlet 
coupling constant <;* 




Figure 3: f y ^'-direction critical temperature vs magnetic field 
of the f y symmetry for t a — 3000K and t c — 5K. T ct ~ 0.5k is 
the critical temperature corresponding to the triplet coupling 
constant gl 



For our numerical calculations we have used the stan- 
dard in-plane transfer integral values, t a : t b ~ 3000 : 
300if [l| . After maximization with respect to the LOFF 
wave vector k, we find H FOFF (t c = 0) = 5.3 T for a 
LOFF wave vector q x = ^^(1 - 0). 

The expected effect of a finite t c is to further reduce 
the value of H^ OFF . In the calculation of the critical 
fields, the values used in the litterature [l^, [H, |H, 
for t c range from 2 K to 17 K, and those used for the 
anisotropy ratio t a /t c range from ~ 200 to ~ 700. How- 
ever, the maximum LOFF critical field, H^ OFF (T = 0), 
depends strongly on the value of the t a /t c ratio Fig. 2] 
Therefore, in the subsequent numerical calculations, we 
will limit ourselves to values for t c in the range 5 — 10K, 
which are the closest values to experimental and theoret- 
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Figure 4: zero temperature critical field dependance on t c . 
The t c ~ 3.8A' value, at which there is a sudden drop of 
H c LOFF (0), corresponds to a semiclassical amplitude A 2 reach- 
ing exactly the interplane distance c, which corresponds to the 
deconfinement of electrons out of the a — b planes. 

ical determinations for t c [2^, Hi| . In that case H^ OFF 
at T = is reduced to the very low value of 0.34T. Then, 
a LOFF state could not be stable and the system is in a 
triplet state at high fields, because of the reentrance phe- 
nomena. However, H^ OFF is so sensitive to t c , which is 
not precisely known, that t c < 3.5JC would be enough to 
get H^ OFF ~ 2.5T. We cannot exclude this possibility, 
which would ensure the existence of a LOFF state and of 
a triplet state at high fields. 

By numerically solving the gap equations pip and 
([12]) , we have been able to calculate T C (H) for the ho- 
mogeneous and LOFF d x 2_ y 2 phases as well as for the f y 
phase Fig. [5] an [3J The most striking result is that, for 
the standard anisotropy ratio t a /t c ~ 600, the maximum 
critical field of the LOFF phase is strongly reduced by 
the anisotropic orbital effects, as can be seen on Fig. [5] 
Therefore, the value of the anisotropy ratio is crucial to 
the determination of the critical field of the LOFF phase. 
In order to reproduce the critical fields measured experi- 
mentally we have to use very high anisotropy ratios which 
do not correspond to the commonly known values. For 
the standard values indeed, < 700, the LOFF phase is 
limited to very low magnetic fields, Fig. Q] This could 
qualitatively explain the recent NMR measurements by 



Shinagawa et al., where at a field of 0.9T the Knight-shift 
in the CIO4 is consistent with singlet symmetry, while at 
slightly higher fields, ~ 1.4T, the Knight shift is that of a 
triplet symmetry. Then, our results support the hypoth- 
esis of a phase transition from singlet to triplet supercon- 
ductivity that takes place at relatively low fields, as can 
be seen on Fig. [5] It should be noted that a first order 
phase transition is also expected from homogeneous to 
LOFF superconductivity [ll|, [T^]. However, the proxim- 
ity of the two transitions, Fig. [5] could make it difficult 
to discriminate between them. 

For the / y -triplet phase, we find the same reentrance 
phenomenon, Fig . E fl as that previously found for the 
p y symmetry [ill IT2T [331 ] . At high magnetic fields, the 
triplet pairing of the Cooper electrons and the absence 
of Zeeman effect are at the origin of the re-stabilization 
of superconductivity. This phenomenon has already been 
observed in the CIO 4 compound by Lee et al. 0. Besides 
the reentrance and the oscillations of T c at low temper- 
atures, our calculations show that, at intermediate mag- 
netic fields, the critical temperature is strongly reduced, 
Fig. [3] This is clearly an orbital effect, and may be 
improved if we take into account the magnetic field de- 
pendance of the coupling constants, which may not be 
negligible in our case. 

Conclusion We have used an RG calculation to de- 
termine the singlet and triplet coupling constants which 
we have used in the calculation of the critical fields. We 
show that for a not too high anisotropy ratio t a /t c , the 
critical fields for a d-wave LOFF state are strongly re- 
duced, while a reentrant superconductivity is expected 
for the f-wave triplet state. From these results we con- 
clude that a transition from the singlet LOFF state to 
the triplet state at low magnetic fields seems to be the 
most reliable interpretation of the recent experimental 
results. 

More experimental work is still needed though. In par- 
ticular, more measurements of the Knight-shift as a func- 
tion of magnetic field at low temperatures are necessary. 
Specific heat measurements will also be needed in order 
to observe the phase transition. 

We would like to thank J. Friedel, J. C. Nickel, S. 
Charfi-Kaddour and S. Haddad for stimulating discus- 
sions and valuable help. 
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Figure 5: Schematic H — T phase diagram calculated from 
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